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ABSTRACT 

The next generation of large sky photometric surveys will finally be able to use 
arc statistics as a cosmological probe. Here we present the first of a series of papers 
on this topic. In particular, we study how arc counts are sensitive to the variation of 
two cosmological parameters: the (total) matter density parameter, ff m , and the nor¬ 
malisation of the primordial power spectrum, expressed in terms of ag. Both these pa¬ 
rameters influence the abundances of collapsed structures and their internal structure. 
We compute the expected number of gravitational arcs with various length-to-width 
ratios in mock light cones, by varying these cosmological parameters in the ranges 
0.1 ^ fi m < 0.5 and 0.6 < as < 1. We find that the arc counts dependence on 
and o 'g is similar, but not identical, to that of the halo counts. We investigate how 
the precision of the constraints on the cosmological parameters based on arc counts 
depends on the survey area. We find that the constraining power of arc statistics de¬ 
grades critically only for surveys covering an area smaller than 10% of the whole sky. 
Finally, we consider the case in which the search for arcs is done only in frames where 
galaxy clusters have been previously identified. Adopting the selection function for 
galaxy clusters expected to be detected from photometric data in future wide surveys, 
we find that less than 10% of the arcs will be missed, with only a small degradation 
of the corresponding cosmological constraints. 
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1 INTRODUCTION 

The importance of galaxy clusters in cosmology is well 
known (for a review, see Allen et al. 2011, and references 
therein). Being the most massive bound systems in the uni¬ 
verse, they trace the latest stage of structure formation. 
Their abundance and mass as a function of redshift are thus 
highly indicative of how the growth of the cosmic density 
fluctuations occurs and can thus be used to constrain the 
matter content, the initial power spectrum normalisation 
and the expansion history of the universe (Eke et al. 1998; 
Borgani et al. 2001; Reiprich & Bohringer 2002; Allen et al. 
2003; Schuecker et al. 2003; Henry 2004; Vikhlinin et al. 
2009; Mantz et al. 2010; Rozo et al. 2010; Sehgal & et al. 
2011; Planck Collaboration et al. 2014a; Benson et al. 2013; 
Mantz et al. 2014). 

Moreover, galaxy clusters are essential cosmic labora¬ 


tories where the complex interaction between baryons and 
dark matter can be studied in detail, in particular during 
merger events (Clowe et al. 2006; Merten et al. 2011). 

Being the most massive structures in the Universe, 
galaxy clusters are also the most powerful gravitational 
lenses (Narayan & Bartelmann 1999; Bartelmann & Schnei¬ 
der 2001; Kneib & Natarajan 2011). In particular, they are 
responsible for highly non-linear lensing effects taking place 
in their densest regions, i.e. in their cores. In this so-called 
“strong” lensing regime, the images of background galax¬ 
ies are heavily distorted, often leading to the appearance 
of gravitational arcs with large length-to-width ratios (see 
Kneib & Natarajan 2011; Meneghetti et al. 2013, and refer¬ 
ences therein). 

The efficiency of galaxy clusters to produce arcs with a 
given ratio between their length l and width w is quantified 
by means of their strong lensing cross section ai/ w . This is 
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defined as the area on the source plane where the source has 
to be located in order to form an arc with such a l/w ratio. 
To have a large strong lensing cross section, the projected 
mass distribution of the lens must be exceptionally dense 
on the plane of the sky, also called “lens plane”. Sometimes, 
this can be the result of projection effects (halos elongated 
along the line-of-sight, superposition of multiple structures 
at different redshifts, etc), but, generally speaking, gravita¬ 
tional arcs trace the highest peaks in the cosmological den¬ 
sity field. Bartelmann et al. (1998) first pointed out that 
counting gravitational arcs may be a competitive method 
to constrain cosmological parameters (see e.g. Fedeli et al. 
2006; Meneghetti et al. 2013). 

Unfortunately, strong lensing events such as gravita¬ 
tional arcs are rare. Given the relatively low number of arcs 
discovered so far and the high inhomogeneity of cluster opti¬ 
cal surveys, the application of arc statistics in cosmology has 
been attempted mainly with the goal of possibly falsifying 
the standard cosmological model (the so-called concordance 
ACDM) rather than actually constraining cosmological pa¬ 
rameters. For about 15 years, scientists have debated on the 
existence or nonexistence of a tension between the observed 
number of arcs and the predictions derived in the framework 
of the cosmological model favoured by observational data 
(see Meneghetti et al. 2013, and references therein). Un¬ 
doubtedly, arc statistics suffered so far from the lack of suit¬ 
able large observational datasets for a reliable comparison to 
theoretical predictions. The situation is likely to change rad¬ 
ically in the near future, thanks to the advent of large optical 
surveys, covering areas in the range from several thousands 
of square degrees (de Jong et al. 2013; The Dark Energy 
Survey Collaboration 2005), to (almost) the full sky (LSST 
Science Collaboration et al. 2009; Karoji 2009; Laureijs et al. 
2011; Spergel et al. 2013). 

In Boldrin et al. (2012) we forecasted the number of 
gravitational arcs visible in the future wide field survey to 
be performed by the ESA Euclid mission (Laureijs et al. 
2011). A further step is the analysis of the sensitivity of arc 
statistics on cosmological parameters. In particular in this 
work, we focus on the dependence of arc statistics on the (to¬ 
tal) matter density parameter f l m and on the normalisation 
of the primordial power spectrum, expressed in terms of ag. 
More precisely, we investigate the region of the parameter 
space defined by = [0.1 — 0.5], and ag = [0.6 — 1.0]. 
We sample the parameters at intervals AQ m = 0.1 and 
Act 8 = 0.1, thus investigating a total of 25 different cos¬ 
mological models, always making the assumption of a flat 
Universe (fi m + S2a = 1). To help a direct comparison with 
the results of Boldrin et al. (2012), here we assume a refer¬ 
ence cosmology defined by the parameters (fl m = 0.272, 
os = 0.809) with present Hubble parameter Ho = 70.4 
km/s/Mpc, in agreement with the WMAP7 results (Ko¬ 
matsu et al. 2011). However we will discuss also the results 
obtained with the assumption of the parameters recently 
derived by Planck Collaboration et al. (2014b). 

The paper is organized as follows. In Section 2 we briefly 
describe the method adopted to compute the number of arcs 
and to build the mock catalogs. In Section 3 we present our 
main results, i.e. the total number of arcs and their redshift 
distribution as a function of Q rn and ag. Section 4 is devoted 
to a discussion of the origin of the cosmological influence on 
arc statistics; we also consider how the results can change 


when including the effects of possible systematics. In Section 
5 we discuss the agreement between the predictions obtained 
with our methodology and the arc counts in real surveys, 
using the recent results of Xu et al. (2015) based on the 
CLASH survey. Finally, in Section 6 we summarise and draw 
our conclusions. 


2 METHOD 

In this section, we summarise the steps undertaken for cal¬ 
culating the number of arcs produced by a population of 
strong lenses. 

2.1 Generation of lenses with M0KA 

In the attempt of properly modelling all the relevant struc¬ 
tural properties of the lenses in our calculations of a t / w , 
we make use of the pseudo-analytic code M0KA (Giocoli 
et al. 2012a, b). This code allows to generate mock cluster- 
size gravitational lenses in any desired cold dark matter 
(CDM) scenario, including features like: triaxiality and pro¬ 
jection effects, scatter in concentration, substructures, and 
the presence of a brightest central galaxy, including the ef¬ 
fects that its growth produces adiabatic contraction on the 
cluster dark-matter halo profile. As discussed in several ear¬ 
lier works based on fully numerical simulations (see e.g. 
Meneghetti et al. 2003c; Torri et al. 2004; Meneghetti et al. 
2007), all these features play an important role for deter¬ 
mining the cluster ability to produce giant arcs. The most 
important novelty of our work mainly lies on the adoption 
of very realistic simulated strong lenses, as the M0KA code 
allows to do. 

More in detail, to generate the mock lenses we use the 
following prescriptions: 

• the host halo mass density profile is described by the 
Navarro, Frenk and White (NFW) radial function (Navarro 
et al. 1996); 

• the dark matter halo concentration is derived from the 
mass-concentration relation proposed by Zhao et al. (2003), 
assuming log-normal scatter of <ti 11c = 0.25; 

• the lens halo is triaxial and the axis ratios are generated 
assuming the triaxial distributions derived by Jing & Suto 
( 2002 ); 

• the spatial orientation is randomly chosen; 

• the substructure abundance follows the sub-halo mass 
function derived by Giocoli et al. (2008); 

• the substructure radial distribution follows that pro¬ 
posed by Gao et al. (2004); 

• the substructures are modelled as Single Isothermal 
Spheres (SISs) (Metcalf & Madau 2001); 

• the mass density profile of the Brightest Cluster Galaxy 
(BCG) resembles the Hernquist’s profile (Hernquist 1990). 

The main outcome of the M0KA code is the deflection an¬ 
gle map of each mock halo. From the latter, we obtain au w 
for sources at z 3 = 2 via ray-tracing technique. For simplic¬ 
ity, we consider only isolated lenses, avoiding the cases of 
merging haloes, even though some studies revealed the im¬ 
portance these events have on the amplitude of the strong 
lensing signal (see e.g Torri et al. 2004). Furthermore, an¬ 
other source of enhancement of the signal is given by the 
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Figure 1. Number of arcs (normalised to the reference WMAP7 cosmology) as function of f2 m and for different values of erg. Different 
panels refer to different redshift bins between z = 0.21 and z = 1.26, as labeled. The reported counts represent the median of 128 different 
light-cone realisations for each combination of the cosmological parameters. Solid blue, dotted cyan, dashed green, long-dashed brown 
and dot-dashed dark orange lines indicate the results for erg = 0.6, 0.7, 0.8, 0.9 and 1, respectively. The results refer to arcs with l/w > 10 
and sources ltr above the mean background noise level. 


large scale structure present along the line of sight, which 
boosts the lens projected mass on the source plane. Since 
beyond purpose of this paper we postpone detailed studies 
of those two aspects to a future work. 

2.2 Number of giant arcs 

We trace light-rays through a uniform grid on the lens 
plane towards the source plane, accounting for the position- 
dependent deflections. The source plane is populated with 
elliptical sources distributed on a regular grid. Their lensed 
shapes are recovered by collecting the light rays hitting 
them. The lengths and widths of the resulting gravitational 
arcs are measured using the method described in several ear¬ 
lier papers (see e.g. Meneghetti et al. 2003c). As in Boldrin 
et al. (2012), the source population has random ellipticity 
and an apparent size depending on the redshift (see Boldrin 
et al. 2012, Fig. 3.). To properly sample the region of the 
source plane, where sources are lensed as giant arcs, we it¬ 
eratively increase the resolution of source grid near the lens 
caustics. Each source is then representative of an area on 
the source plane defined by the local resolution of the source 
grid. Using this area as weight, we compute au w , defined as 
the area on the source plane within which a source has to 
be located for being lensed as a giant arc with a given l/w 
(Meneghetti et al. 2008). 


The number of arcs with length-to-width ratio larger 
than l/w, produced by a lens of mass M at redshift zi, can 
be computed solving the following relation: 

r z s,max 

N i/w (M,zi) = a l/w {2) dz s fa{zi,z s )n(zs,S ), (1) 

J z t 

where ai/ w ( 2) = zi, z 3 = 2) is the strong lensing 

cross section for giant arcs for sources at redshift z s = 2, and 

f a = - . accounts for the scaling of oii w with z s . 

&l/w \%s = 

In order to minimize the computational time, we recover the 
scaling function f a computing ai/ w (z 3 ) of a lens subsample, 
adopting different values of z s (see for more details Boldrin 
et al. 2012). 

A crucial ingredient is the redshift distribution of the 
sources exceeding a given surface brightness S, n(z 3 , S). We 
derive it by simulating an observation of the galaxies in the 
Hubble-Ultra-Deep-Field (HUDF) using the SkyLens code 
(Meneghetti et al. 2008), to a depth which is reasonable 
for a future wide field survey from space. As a reference, 
we consider the Euclid wide survey, which is expected to 
reach an average limiting magnitude m r ; z = 24.5 (Laureijs 
et al. 2011). Adopting typical Euclid-like exposure times and 
background levels, we estimate the surface brightness limits 
corresponding to ler and 3<r above the mean background 
level, and use the simulated observations of the HUDF to 
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Figure 2. Number density of arcs as a function of for different values of as- The right and left columns refer to sources detectable 
at lcr and 3 a above the mean background noise level, respectively. From top to bottom, the different panels show the results for three 
choices of minimum l/w, namely 5, 7.5 and 10. Line and color styles are as in Fig. 1. In each panel the horizontal dashed line shows the 
counts in the considered reference mode.. 
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Figure 3. The abundances of arcs with l/w ^ 5 (left panel) and l/w ^ 7.5 (right panel) relative to the abundances of arcs with l/w ^ 10 
as a function of Q m . Line and color styles are as in Fig. 1. The results refer to arcs detectable at lcr above the mean background level. 


measure n(z s , S). Using a subset of MOKA lenses to simulate 
observations of gravitational arcs, we find that we can safely 
neglect from the giant arc counts sources above z 3t max = 6. 

For each combination of cosmological parameters, we 
produce a catalog of cluster-sized lenses with different 
masses and redshifts. We define 100 mass bins which are uni¬ 
formly spaced in logarithm between 10 13 and 10 16 
and 8 redshift bins, having Az = 0.03 and centered at red- 
shifts 0.21, 0.36, 0.54, 0.6, 0.84, 1.02, 1.14, and 1.26. The 
choice of such redshift bins is optimised for the expected red- 
shift distribution of the lenses producing giant arcs, which 
we derived in Boldrin et al. (2012) for a Euclid-like survey. 
For each combination of redshift and mass we use MOKA to 
generate 100 halos with different structural properties and 
measure their ai/ w , from which we can derive the number 
of giant arcs they produce, as discussed above. 

The catalog of lenses is then used to generate 128 re¬ 
alisations of lens distributions (light-cones). In each light- 
cone, which subtends an area of 15,000 square degrees, we 
calculate the number of lenses of mass M and redshift zi ac¬ 
cording to the Sheth & Tormen (1999) mass function, and 
estimate the total number of arcs by summing up the contri¬ 
butions from each individual lens. Finally, we combine the 
different light-cones to measure the median number of arcs 
per square degree and the relative scatter as a function of 
the considered cosmological parameters. 

3 RESULTS 

3.1 Number of arcs as a function of redshift 

We begin by discussing how the arc counts change as a func¬ 
tion of cosmology in different redshift bins. In Fig. 1 we show 
the number of arcs, normalized to the reference WMAP7 
cosmology, as function of Q m . The different panels refer to 
the eight redshifts where the calculations were performed. 
Different colors and line styles are used to display the re¬ 


sults for several values of a g: solid blue, dotted cyan, dashed 
green, long-dashed brown and dot-dashed dark orange lines 
refer to crs=0.6, 0.7, 0.8, 0.9 and 1, respectively. Long-dashed 
black horizontal lines correspond to unity, i.e. to the refer¬ 
ence cosmology. The lack of a blue solid line in the last panel 
is due to the inefficiency of clusters at zi = 1.26 to produce 
giant arcs in the cosmology with as = 0.6. As expected, at 
all redshifts, the arc counts grow as a function of S2 m and 
as a function of as, indicating that the abundance of giant 
gravitational arcs is higher in cosmological models with more 
matter and higher normalisation of the power spectrum of 
the primordial density fluctuations. 

We also notice that the change of arc counts as a func¬ 
tion of cosmology depends on the lens redshift. The depen¬ 
dence on is stronger at lower redshift, and flattens off 
as zi increases. On the contrary, it appears that the value of 
a 8 affects the results more significantly at high redshift. 

While the results in Fig. 1 refer to arcs with l/w ^ 
10 and sources above the lcr background level, the trends 
remain similar for other l/w ratios and detection limits. 


3.2 The total number of arcs 

From the distributions obtained from the 128 different light- 
cone realizations, we measure the median number of arcs per 
square degree expected in each cosmological model. This has 
been done by performing a spline interpolation through the 
above-mentioned 8 redshifts up to a maximum lens redshift 
of zi = 1.5. 

In the reference WMAP7 cosmology, the expected num¬ 
ber densities of arcs per square degree with l/w ^ 5, 7.5, and 
10 are 0.594±0.016, 0.194±0.006, and 0.085±0.003, respec¬ 
tively. These are in excellent agreement with our estimates 
reported in Boldrin et al. (2012), although these were ob¬ 
tained using a larger number of redshift bins and avoiding 
the interpolation. 

In Fig. 2, we show the median arc number counts per 
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square degree as a function of O m . We also show how the 
counts vary by changing the value of erg, using the same color 
and line styles used in Fig. 1. From top to bottom, we show 
the results for l/w ^ 5, 7.5, and 10, respectively. The left 
and the right panels refer to detections at 1- and 3a above 
the level of the background. Obviously, the results show the 
same dependence on and as reported in Fig. 1. 

We also see that the ratios between counts of arcs with 
different l/w depend on the cosmological parameters. As 
shown in Fig. 3, for low I2 m , the abundance of arcs with 
l/w ^ 5 or l/w ^ 7.5, relative to that of arcs with l/w ^ 10, 
is higher, indicating that halos in these cosmological mod¬ 
els have smaller critical lines and are thus less efficient at 
producing large distortions. The ratios also depend on as; 
in cosmologies with higher as halos are able to produce a 
higher abundance of arcs with large l/w. The results (here 
shown only for sources lcr above the mean background level) 
are insensitive to the assumed detection limit. Therefore, in 
the following discussion we will show the results only for arcs 
detectable at the lcr level. We will also focus on arcs with 
l/w ^ 10. 

The upper panel in Fig. 4 shows the difference in the 
arc counts relative to the reference WMAP7 cosmology in 
the — as plane. Within the ranges explored in this work, 
we may find differences of up to one order of magnitude 
for the predicted arc counts between cosmological models. 
We also notice that the cosmological parameters and 
as are degenerate with respect to the arc counts. Indeed, 
the same number of arcs is expected in cosmologies whose 
combination of and erg lays in a banana-like region ex¬ 
tending from the upper left to the bottom right corner of the 
plane. The origin of this degeneracy will be better discussed 
in Section 4.1. Interestingly, a Planck-like cosmology with 

= 0.3086 and erg = 0.8288 (Planck Collaboration et al. 
2014b) produces 54% more arcs than the reference WMAP7 
cosmology. 

We find that in the case of the reference WMAP7 model, 
the equation describing the degeneracy curve between the 
cosmological parameters has the following form: 


O°0.8 




— Aas -f- Bag -(- C , 


( 2 ) 


where A = 1.771, B = —3.952 and C = 2.31. Such function 
is given by the white line in the upper panel of Fig. 4. 

In the attempt to quantify the uncertainty in the arc 
counts, we define the la uncertainty on the number counts 
as a = (a^y + a%) 1 ^ 2 , where acv is the cosmic variance, 
which is estimated from the 16th and 84th percentiles of the 
distributions derived from the 128 light-cone realisations of 
each tested cosmological model. The other term appearing 
in the equation, ap = VdV, is the associated Poisson noise 
on the number counts. 

In the bottom panel of Fig. 4, we perform an error anal¬ 
ysis showing the levels corresponding to 1, 3, and 5 a devia¬ 
tions (from dark to light colors) from the WMAP7 and the 
Planck cosmologies in the — erg plane. The results were 
obtained assuming a survey covering 15,000 sq. degrees of 
the sky to the depth expected for the Euclid mission. It is 
interesting to notice that a survey with the Euclid charac¬ 
teristics will be able to distinguish these two cosmological 
models at more than 5<r level. 


Figure 4. Upper panel: difference in the arc counts with respect 
to the reference WMAP7 cosmology in the — erg plane. The 
results are shown for arcs with l/w s 10 detectable at lcr above 
the background level. The white crosses represent the cosmologi¬ 
cal models having the same arc counts as the reference WMAP7 
model [relation (2)]. Bottom panel: levels corresponding to 1, 3, 
and 5(7 deviations (from dark to light colors) from the WMAP7 
(blue) and the Planck (yellow) cosmologies in the — erg plane, 
assuming a 15,000 deg 2 survey to the expected depth of the Euclid 
wide survey. The crosses indicate the position of the two reference 
models. 


4 DISCUSSION 

4.1 Influence of the cosmological parameters on 
arc statistics 

In this section we will discuss in more detail some aspects 
of the influence of and a g on arc statistics. In general, 
the cosmological parameters play an important role in arc 
statistics through the lens mass function and their strong 
lensing cross section, the latter depending on the geometry 
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Figure 5. The strong lensing selection function (black solid 
curve), i.e. the minimum galaxy cluster mass expected to pro¬ 
duce critical lines for sources located at z s = 2 (Meneghetti et al. 
2010a; Boldrin et al. 2012). For comparison, the red dashed curve 
represents the minimum mass of galaxy clusters which are ex¬ 
pected to be detected above three times the rms of the field galaxy 
counts in the Euclid photometric survey (?). 
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Figure 6. Number density of expected strong lenses as a function 
of redshift, for cosmologies with different Q m and as- Plots from 
left to right and from up to bottom refer to increasing values of 
Qm- Different colors represent counts for various values of as, as 
labeled on the bottom right. The black line shown in all panels 
represents the results for the reference WMAP7 model. 


Halo Counts 



Figure 7. Relative difference of halo counts on the crg-fi m plane 
with respect to the reference WMAP7 model. The white crosses 
represent the degeneration curve relative to the arc counts, while 
the white solid line refers to the degeneracy curve for halo counts. 

of the Universe and on the structural properties of the lens 
halo. 

In particular, the number of arcs is directly related 
to the number of lenses able to produce arcs. Following 
Meneghetti et al. (2010a, 2011), this can be estimated in¬ 
cluding in the mass function describing the lens distribu¬ 
tion a sharp cut at the minimum mass corresponding to the 
smallest systems in which we expect to find critical lines for 
sources at z s = 2. The shape of the adopted selection func¬ 
tion as a function of redshift is shown by the black curve in 
Fig. 5 (see also Boldrin et al. 2012). 

In Fig. 6, we present the number density (given per 
square degree) of the lenses as a function of redshift. In 
each panel, we keep fixed S2 m as labeled and we vary the 
value for as, using the color code indicated on the bottom 
right. To facilitate the comparison, the lens number density 
in the reference WMAP7 cosmology is shown in black in 
all panels. From the figure, the strong effect of the different 
matter density on the lens abundances and the anticipated 
structure formation originated by a higher power spectrum 
normalization are clear. 

In Fig. 7, adopting the same color code as in the upper 
panel of Fig. 4, we show the difference in the lens counts 
relative to the reference WMAP7 cosmology in the — erg 
plane. The white solid curve in the figure represents the 
degeneracy between and as for the halo counts, for which 
we find the following relation: 

<Tg(S2 m /0.272) 0 ' 304 = 0.809 . (3) 

Even if with some differences, this curve is close to the re¬ 
lation (shown by the white crosses) representing the degen¬ 
eracy we found in the f l m — as plane for the arc counts (see 
also Fig. 4): this is clearly due to the fact that the most 
important ingredient for arc statistics is the lens mass func- 
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lens redshift 

Figure 8. Lensing distance for different values of f! m . Sources 
are kept fixed at redshift z s = 2. 

tion. However, if one compares the amplitude of the count 
variation by looking at the width of the coloured strips, it 
is evident that the arc density is more sensitive to the cos¬ 
mological parameters than the simple halo density: conse¬ 
quently, a wide survey of gravitational arcs could potentially 
give significant constraints. 

The larger sensitivity of arc statistics is due to the cos¬ 
mological dependence of the other main ingredients, such as 
the angular diameter distances of lenses and sources and the 
lens structural properties. We know that the first condition 
for an axially symmetric lens to act like a strong lens is that 
in some points x on the lens plane the condition 

k(x ) > 1 (4) 


at redshift z s = 2 and we study Di en3 (zi), that is we keep 
fixed the length of the lensing system and we move the lens 
from the observer towards the source plane. The results are 
shown in Fig. 8: we see that increasing the value of the 
strong lensing efficiency reaches its maximum at lower lens 
redshifts. In particular the peak around which the produc¬ 
tion of gravitational arcs is expected to be boosted shifts 
from z « 0.6 to z ~ 0.4 when the value of is increased 
from 0.1 to 0.5. 

The effect of the anticipation of structure formation due 
to a higher value of as (Giocoli et al. 2007, 2012b) has con¬ 
sequences on several halo structural properties that may in¬ 
fluence the size of ai/ w . Considering the concentration pa¬ 
rameter, at fixed as, large fi m values lead to larger con¬ 
centrations because the structures form and grow in denser 
environments. At the same time, keeping fixed the value of 
flm, in cosmologies with high as the concentration increases 
because of both the higher contrast between primordial per¬ 
turbations and background, and the anticipated formation 
time (Neto et al. 2007; Giocoli et al. 2012b; Maccio et al. 
2008). 

Haloes triaxiality is also an important feature that 
depends on cosmological parameters (Despali et al. 2014; 
Bonamigo et al. 2015). In particular, the level of spheric¬ 
ity of a halo, which is directly related to the ratio between 
its minor and major semi-axes o/c, is an increasing func¬ 
tion of ag and a decreasing function of fi m . As an exam¬ 
ple, if we consider haloes with a mass equal to 7.5 x 10 14 
ft- 4 M 0 at redshift z = 0.54 in a cosmological model with 
Qm = 0.3, the median ratio among 128 realizations varies 
from a/c = 0.353lQ Qgg in a model with erg = 0.6, to 
a/c = 0.417lQ Q|g in a model with <T8 = 1.0. The quoted 
uncertaintes correspond to ler errors. On the other hand, if 
we fix as = 0.8, the ratio changes from a/c = 0.419lg'Qgg 
in a model with Q m = 0.1, to a/c = 0.3881q'q|i in a model 
with = 0.5. 


4.2 Effects of completeness and cluster selection 
function 


occurs, where k = E(a?)/E cr is the so-called convergence, 
E(x) is the lens projected mass density and 


E cr 


4ttG 


D 


-i 

lens 


(5) 


represents the critical value of the two-dimensional mass 
density in order to have strong lensing effects. The quan¬ 
tity Di e n S is the so-called lensing distance, defined as 


D 


lens — 


DlsDl 

D s 


( 6 ) 


where Ds, Dl and Dls are the angular diameter distances 
of the source, of the lens and between source and lens, re¬ 
spectively. Although for elliptical lenses we have to add the 
effect of shear to the condition (4), we can infer, to first 
approximation, what are the system configurations which 
are more efficient in producing strong lensing features by 
investigating how Di ens changes in the different cosmologi¬ 
cal models, once the lens properties and the source redshifts 
are kept fixed. We remind the reader that Di ens contains the 
full dependence on the geometry of the system and does not 
depend on os, but only on S2 m . We fix the source position 


In the following subsections we will discuss how our results 
change when we take into account the lack of completeness 
and when we introduce a realistic photometric galaxy cluster 
selection function. 


f.2.1 The effect of sample completeness and survey area 

Let us consider here the case in which a fraction of arcs are 
missed, independently of the properties of the lens configura¬ 
tions (l/w, zi, z 3 ). This may happen because some arcs may 
escape detection for some particular configurations of the 
light distribution within the cluster, or when the separation 
between cluster and foreground galaxies is made difficult by 
the lack of precise color information. The total arc counts 
may also diminish because we are performing our search in 
a reduced effective area, smaller than the one of the run¬ 
ning survey. In this situation losing 10% of the counts is 
equivalent to observe a portion of sky 10% smaller than the 
original survey. The obvious consequence of a reduction of 
the number of observed arcs is that the Poissonian uncer¬ 
tainty grows and can start to dominate with respect to the 
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Figure 9. The amplitude of the 3 a error bar on f2 m as a function 
of the survey area. The vertical dashed line shows the size of the 
future Euclid wide survey. The red horizontal dotted and dashed 
lines represent a variation of corresponding to ±10% and 
±25%, respectively. 


Figure 10. Percentage of arcs effectively detected by considering 
only lenses having a mass larger than the Euclid cluster photo¬ 
metric selection function. Different color refer to different values 
of as , as labeled; solid and dotted lines are for arcs with l/w ^ 5 
and 10. 


cosmic variance, when accounting for the total error budget. 
To quantify this effect, in Fig. 9 we show, as a function of the 
fraction of the sky covered by the arc search, the variation 
of the 3cr error bar on the parameter fi m , when the value 
of as is a priori fixed to its reference value (os = 0.809), as 
it may happen if independently measured from other cos¬ 
mological probes. Dark, medium and light blue regions refer 
to the cases of arcs with l/w ^ 5, 7.5 and 10, respectively, 
while the horizontal dotted (dashed) lines indicate an accu¬ 
racy of 10 (25) per cent on O m . From the figure it is clear 
that arcs with l/w ^ 5, being more numerous, give stronger 
contraints and are less affected by possible incompleteness 
problems. However, there is a difficulty when dealing with 
them because they can look like simple edge-on galaxies. 
For this reason the loss and misidentification of arcs are 
expected to depend on l/w, being stronger for low -l/w ra¬ 
tios. From this point of view, Fig. 9 is quite encouraging: if 
the survey area is sufficiently wide (larger than 10% of the 
whole sky), or equivalently if the arc finders are sufficiently 
efficient, the error budget is dominated by cosmic variance 
and there is not a significant difference in the constraining 
power between using arcs with l/w ^ 5 or with l/w ^ 10. 
We remind that the SDSS (York et al. 2000) has an area of 
about 10,000 deg 2 , while the Euclid wide survey is expected 
to cover 15,000 deg 2 (Laureijs et al. 2011). 


4-2.2 The effect of the cluster selection function 

Due to the high computational cost of the algorithms for arc 
detection, a possible strategy in future wide surveys is to run 
these codes only on small-size frames where galaxy clusters 
have been previously identified. Obviously, this originates a 
reduction of the effective number of arcs, which is strongly 



Figure 11. As the upper plot of Fig. 4, but for arcs produced by 
lenses having a mass larger than the Euclid cluster photometric 
selection function (?). The crossed line represents the degeneracy 
curve obtained when no selection function is applied. 

dependent on the specific cluster selection function of the 
survey. 

As a worked example, here we consider again the future 
ESA Euclid mission. Given the amount and quality of its 
data covering an area of 15,000 deg 2 , there will be at least 
three main ways to identify galaxy clusters: (i) from pho¬ 
tometric data, (ii) from spectroscopic data, and (iii) from 
cosmic shear maps. As shown in ? , the one based on photo- 
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metric data (see, for example, Bellagamba et al. 2011, and 
references therein) is expected to be largely the most efficient 
one. In this case, the minimum mass of galaxy clusters hav¬ 
ing a number of members larger than 3 times the r.m.s. of the 
field galaxy counts is expected to be between 5 x 10 13 Mq/ h 
and 8 x 10 13 Mq//i in the redshift range here considered (?). 
Compared to the minimum mass needed to produce critical 
lines for sources located at redshift z s = 2 (see Fig. 5), the 
Euclid cluster selection is then slighly higher on a limited 
redshift range only, namely between z = 0.2 and 2 = 0.5. 
This means that limiting the search for arcs to frames where 
galaxy clusters have been already identified is expected to 
not reduce dramatically the number of detected arcs. This is 
confirmed in Fig. 10, where we show the fraction of arcs that 
can be effectively detected following this strategy. Same col¬ 
ors indicate same values of as, as labeled in the figure, while 
solid and dotted lines refer to arcs with l/w ^ 5 and 10, re¬ 
spectively. For the cosmological models here considered, the 
reduction varies between 2 and 10 per cent and is almost 
independent of l/w. For the reference WMAP7 model, the 
percentage of effectively detected arcs remains about 95 per 
cent. 

In Fig. 11 we show the relative differences in the arc 
counts between each cosmological model and the refer¬ 
ence WMAP7 cosmology, considering only arcs produced 
by galaxy clusters above the Euclid photometric selection 
function. In the figure, the color scale is identical to that 
adopted in the upper panel of Fig. 4. The white crosses 
represent the degeneration curve we found considering the 
total number of arcs, i.e. without applying the cluster se¬ 
lection function. Although similar, the curve changes in a 
non negligible way, especially considering extreme values of 
the parameters. This underlines the importance of taking 
into account every kind of selection function when combin¬ 
ing theory and observations in arc statistic studies, avoiding 
possible systematics. 

The presence of strong lensing features like arcs can 
represent a complementary way to confirm the presence of 
a galaxy cluster. Moreover arcs can be used to improve the 
estimates of the mass of galaxy clusters, a fundamental in¬ 
gredient to fully exploit the evolution of their abundance 
as cosmological probe. For this reason it is important to 
compute what is the fraction of the galaxy clusters iden¬ 
tified in the Euclid photometric survey, which are able to 
produce at least one giant arc. The result for the reference 
WMAP7 cosmology as a function of redshift is shown in 
Fig. 12 for arcs with l/w ^ 5, 7.5 and 10 (black solid, red 
dashed and green long-dashed lines, respectively). Typical 
mean values are around 1 per cent, 0.33 per cent and 0.15 
per cent for l/w ^ 5, 7.5 and 10, respectively. From the 
figure we notice that the strong lens fraction peaks around 
2 = 0.5: this behaviour is a combined effect between the 
well of the photometric selection function around redshift 
z = 0.75 and the peak - around the same redshifts - of the 
strong lens counts. Interestingly, for redshift z ^ 1.3 the per¬ 
centage tends to vanish. Finally we notice that the fact that 
the strong lensing selection function can be smaller than the 
photometric cluster one would allow in principle to add ex¬ 
tra objects to the Euclid cluster sample by looking for strong 
lensing features only. However this would require to run the 
algorithms for arc detection blindly in different areas of the 
survey. Considering the reference WMAP7 model and arcs 



lens redshift 

Figure 12. Fraction of galaxy clusters having a mass larger 
than the Euclid cluster photometric selection function produc¬ 
ing at least one giant arc. Results are shown for the reference 
WMAP7 model. Different line styles refer to different length-to- 
width ratios, as labeled. 

with l/w ^ 5, the gain would correspond to approximately 
300 extra objects only, all having a relatively low redshift 
(0.2 ^ 2 ^ 0.5). 

Therefore, we can conclude that arc statistics represents 
a complementary tool to identify galaxy clusters or eventu¬ 
ally to prove their presence. In particular, arcs with a small 
l/w ratio are the best tracers, since they are more numerous, 
but, at the same time, they are the more difficult to identify 
because of their similarity with non-lensed galaxies. Finally, 
our results underline that the codes for arc identification 
can be run on single frames where galaxy clusters have been 
already detected with no consequences on the cosmological 
predictive power of arc statistics. 


5 A TEST-BED FOR THE METHOD: THE 

CLASH SURVEY 

While this paper focuses on the sensitivity of arc statistics to 
cosmological parameters like and erg, it is worth mention¬ 
ing that another paper has been recently submitted by our 
collaborators Xu et al. (2015) to compare theoretical predic¬ 
tions of arc abundances in a ACDM cosmological model and 
observations. More precisely, in this other work MOKA has 
been used to build up halos reproducing the properties of 
the X-ray selected galaxy clusters belonging to the CLASH 
sample (Postman et al. 2012). Numerical hydro-dynamical 
simulations tailored to reproduce the CLASH selection func¬ 
tion (Meneghetti et al. 2014) are also used to derive theoret¬ 
ical predictions. Thus, the work of Xu et al. (2015) provides 
the best opportunity for validating our methodology against 
more complex models of the cluster mass distribution and 
against observed clusters with a known selection function. 

The results of this study show that there is an excellent 
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agreement between expectations based on MOKA halos and 
numerical simulations and the arc counts in the CLASH 
clusters. More specifically, the lensing efficiency measured 
in the CLASH sample is 4 ± 1 arcs (with l > 6” and 
l/w >7) per cluster. MOKA simulations return exactly the 
same number (4± 1), while numerical simulations give 3± 1 
arcs per cluster. Therefore, according to Xu et al. (2015), in 
terms of efficiency to produce long and thin arcs, observa¬ 
tions and simulations based on MOKA and numerical hydro- 
dynamical techniques come into full agreement. It is partic¬ 
ularly significant that the methodology we have developed 
for modeling cluster lenses for arc-statistics calculations is 
fully capturing the complexity of numerically simulated ha¬ 
los, as evinced from the fact that the cross sections for giant 
arcs of MOKA generated halos are well matching those of 
the halos described in Meneghetti et al. (2014). 


6 CONCLUSIONS 

In this work we have investigated how the number of grav¬ 
itational arcs depends on cosmology, focusing our attention 
on the (total) matter density parameter 11 rn and on the ini¬ 
tial normalisation power spectrum parameter as- In more 
detail we have considered the ranges fl m = [0.1 — 0.5] and 
as = [0.6 — 1.0]. Our main results can be summarised as 
follows. 

• We confirm that arc statistics is very sensitive to the 
couple of parameters fi m — a g. In particular we find that the 
expected number of arcs is an increasing function of both 
parameters: this is mostly due to the fact that increasing 
these parameters boosts the number of lenses. 

• The efficiency in producing arcs in cosmologies with 
high values of as is larger, since it has an effect also on 
the structure formation time, that in turn affects some lens 
structural properties (mainly concentration and triaxiality) 
relevant for strong lensing. 

• A strong degeneracy exists between the two considered 
cosmological parameters for the number of arcs Narcs', for 
the reference WMAP7 model this is expressed by the re¬ 
lation (2), that is similar, but not equal, in shape to the 
degeneracy derived from galaxy cluster counts (see eq. 3). 
The differences between the two arise from the non negligible 
contribution to ai/ w given by the lens structural properties 
- triaxiality, asymmetries, concentration, substructures and 
the BCG - and the lensing distance relation. 

• Arcs with small l/w ratio are more suitable to constrain 
cosmological parameters, since they are more numerous. On 
the other hand, they could be more difficultly identified in 
the surveys because of their similarity with non-lensed galax¬ 
ies. We find that if the survey area is sufficiently larger (more 
than 10% of the full sky) the error budget is dominated by 
cosmic variance, and the constraining power of arc counts 
becomes almost independent of the value of l/w. In particu¬ 
lar a survey covering 15,000 deg 2 will be able to distinguish 
at more than 5 ct level the two cosmological models supported 
by WMAP7 and Planck CMB data. 

• Considering future wide surveys, like the ESA Euclid 
mission, we find that searching for arcs only in frames where 
galaxy clusters have been previously detected will produce 
a loss of 2-10% of arcs only (depending on the cosmological 


model) and a consequent limited degradation of the con¬ 
straining power of arc counts. This suggests that it will be 
not necessary to run the computationally expensive algo¬ 
rithms for arc detection on whole wide surveys. 

In this paper we have discussed the potentiality and the 
capability of the giant arc statistic to constrain the matter 
density and the initial power spectrum normalisation param¬ 
eter in light of the large data-sets that will become available 
from future wide-held surveys. 
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